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Abstract 

We show how to use boundary conditions to drive the evolution on 
a Quantum Mechanical system. We will see how this problem can be 
expressed in terms of a time-dependent Schrodinger equation. In par¬ 
ticular we will need the theory of self-adjoint extensions of differential 
operators in manifolds with boundary. An introduction of the latter as 
well as meaningful examples will be given. It is known that different 
boundary conditions can be used to describe different topologies of the 
associated quantum systems. We will use the previous results to study 
how this topology change can be accomplished in a dynamical way. 

Keywords: Boundary dynamics; topology change; domain depen¬ 
dent Schrodinger equation; quantum Faraday law. 


1 Introduction 


A first motivation for the study of quantum systems with boundary is that 
it is almost necessary to deal with the boundaries of the system when one 
is describing the interaction of a system with the outer world. Even if the 
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system is considered to be isolated at some region in space, one necessarily 
needs to deal and describe the interaction of the system with its walls. In 
addition to this, quantum systems with boundary are getting increasing 
attention from the scientific community since several interesting physical 
phenomena are related with the presence of boundaries, like Casimir effect, 
Quantum Hall effect or Topological insulators [1, 2, 3, 4], 

What will be more important for us in the context of the present arti¬ 
cle is that boundary conditions are an effective way of describing different 
topologies. We will devote the rest of this section to give a brief introduction 
to quantum mechanics which will also serve to hx the notation. In Section 2 
we will introduce the Laplace-Beltrami operator and will show that different 
boundary conditions do characterise different self-adjoint extensions of this 
operator. This identification is used in Section 3 to set the ground for a 
proper description of dynamical change of the boundary data. Finally, in 
Section 4 we address the problem of changing the topology of a quantum 
system in a dynamical way. 


The space of states in quantum mechanics is provided by a complex 
separable Hilbert space Ti.. The scalar product is denoted by (•, •) and || • || 
stands for the corresponding norm. Later on we will discuss time-dependent 
problems. To establish the main notions needed to address them, let us first 
introduce the time-independent problems. Dynamics in Quantum Mechanics 
is governed by the Schrodinger equation: 

= H^{t ), 

where H is a linear, self-adjoint, generally unbounded operator acting on the 
Hilbert space of the system. This operator is called the Hamiltonian opera¬ 
tor. If it does not depend explicitly on the time t, Schrodinger’s equation is 
said to be time-independent. Solutions of the time-independent Schrodinger 
equation can be obtained by means of strongly continuous one-parameter 
unitary groups. 

Definition 1. A strongly continuous one-parameter unitary group is a one 
parameter family of unitary operators [7 : M —)• U{T-L) such that 

U{ti)U{t 2 ) = U{ti + t 2 ) , U{0)=ln, 

lim||[/(t)4'-^|| =0, y'i/en. 
t—>-0 

Stone’s theorem establishes a one-to-one correspondence between strongly 
continuous one-parameter unitary groups and self-adjoint operators. This 
correspondence can be summarised by the expression; 

U(t) = exp(itiL) . 
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Given an initial state 'to £ T^{T~L)i a solution of Schrodinger’s equation is 

^(t) = C/(t)^'o. 

The operator H is generally an unbounded operator. On the contrary to 
what happens in the bounded case, there is a difference between symmetric 
and self-adjoint operators. We will illustrate this difference by means of the 
following example. 

Example 2 (The momentum operator on the interval [0,1]). We will con¬ 
sider the differential operator defined by 



acting on the Hilbert space of square integrable functions T^([0,1]). Notice 
that one needs to select a domain for this operator since there are functions 
T G T^([0,1]) that do not verify PT G T^([0,1]). For instance, T = ^/x . A 
suitable choice for the domain is 'P(P) = ((0,1)), the space of smooth 

functions with compact support contained in the interior of the interval. 

Now we want to obtain the domain of the adjoint operator. In order to 
determine it we will use Green’s formula: 

(T , Pd>) - (PT , $) = i [T(l)$(l) - T(0)d>(0)] . 

It is clear that the right hand side vanishes no matter what the values of 
'I'(l) and 'I'(O) are, since G ((0,1)) . Hence, the only requirement for 
T to be in the domain of the adjoint operator is that PT G T^([0,1]), i.e., 
that the first derivative must be also a square integrable function. The space 
of such functions is called the Sobolev space of order one and is denoted as 
P^([0,1]). It is clear now that 

P(pt) = ^i([0,1]) D ((0,1)) = P(P). 

Thus, the momentum operator P with domain ((0,1)) is symmetric 
but not self-adjoint. Being symmetric means that the identity 

(T,P^>) - (PT,$) = 0 

holds for all T G P(P) • To be self-adjoint one needs the extra require¬ 
ment that 'D(Pf) = P(P), what is not satisfied under the conditions of the 
example. That P(P^) D P(P) is not a particular property of the example 
but the general situation for symmetric operators. If one wants to select a 
domain in which the momentum operator is self-adjoint one needs to select a 
larger domain Vs.a. P(P) such that P(P'I’) D Pj.a. = D P(P) • The 
general characterisation of those operators admitting such a construction 
was done by von Neumann, cf. [5] and references therein. For differential 
operators however, one shall take a more direct approach in terms of bound¬ 
ary conditions. This can be proved to be equivalent to the former theory 
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[6]. It is then necessary to select a space that makes the boundary term 
coming from Green’s formula vanish. For instance, one can take the family 
of domains 

V{Pa) = {^en^ [0,1] I $(0) = e*“$(l)} , a G [0, 27r) . 

For each value of the parameter a the domain above defines a different 
self-adjoint extension of the momentum operator. 


2 Self-adjoint extensions of the Laplace-Beltrami 
operator 

We shall consider the dynamics of a quantum non-relativistic free particle 
moving on a Riemannian manifold. In particular we will consider that the 
manifold has boundary, so that different topologies can be described. Let 
0 be a smooth, compact, Riemannian manifold with smooth boundary 50. 
The dynamics of a quantum free particle constrained to the manifold O are 
governed by the Laplace-Beltrami operator. If we denote the Riemannian 
metric by rj, this operator can be defined in terms of local coordinates {x*} 
as 

The Hilbert space of the system shall be the space of square integrable func¬ 
tions on the manifold £^(0). As in the case of the momentum operator we 
can use the boundary data to characterise the different self-adjoint exten¬ 
sions. We use small size greek letters to denote the corresponding restriction 
to the boundary ip := . We use dotted small greek letters to denote 

the restriction to the boundary of the normal derivative pointing outwards 
If := Green’s formula for the Laplace-Beltrami operator can be 

written using this notation as 

($ , -AT) - (-AT , T) = (v? , iI>)qq - {ip , , 

where {■ ,-)dn denotes the scalar product of the Hilbert space of square 
integrable functions at the boundary, . To determine boundary con¬ 

ditions that lead to self-adjoint extensions of the Laplace-Beltrami operator 
we need to look for maximal subspaces that make the boundary term van¬ 
ish. The maximality requirement is needed to ensure that the extension is 
truly self-adjoint and not just symmetric. One can obtain these boundary 
conditions by means of the Gayley transform, cf. [7, 8], defined by 

p+ := p + ip , p- := p — ip . 
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In terms of these new variables the boundary term becomes 

, V’+)ao - i{^- , i’-)dn ■ 

It is clear that boundary conditions of the form 


(p- = U(p+ (p — i(p = U{(p + i(p) , (1) 

where U £ U [C‘^{dVL)) is a unitary operator on the Hilbert space of the 
boundary, make the boundary term vanish. That these are maximal and that 
every maximal vanishing subspace of the boundary form can be expressed in 
this way is a well known result [9, 7, 8 ]. One-to-one correspondence between 
self-adjoint extensions and unitary operators in the Hilbert space of the 
boundary can only be achieved in dimension one. Fortunately, imposing 
extra regularity conditions on the unitary operators, one can still describe 
a wide class of self-adjoint extensions of the Laplace-Beltrami operator. In 
particular, those exhibited in Ex. 5, see [ 8 , 10]. 

Example 3 (Free particle on the interval H = [0,1]). The boundary of the 
manifold is in this case = { 0 , 1 } and therefore the unitary group acting 
on the Hilbert space of the boundary is U ((£^(9H)) =U{2). The simplest 
cases are when U = — 112 x 2 and U = 12 x 2 , that correspond respectively with 
Dirichlet and Neumann boundary conditions. If one considers 

o) ■ ® 

the boundary condition of Eq. (1) becomes 

^.( 0 ) = $(I) , $'( 0 ) = ^'(I) . 

That is, periodic boundary conditions joining both ends of the interval. On 
the other hand if one takes 



the boundary condition of Eq. (1) becomes 

$(0) = e*“^>(l), ^>'(0) = e*“$'(l) . 


(3) 


We call the latter quasi-periodic boundary conditions. 

Example 4 (Free particle in two intervals, Q = Ii U I 2 ). In this case the 
boundary is the set dQ = {ai, 62 , 02 , 62 } (see Fig. I), and the unitary group 
at the boundary is therefore U [C‘^{dQ)) ~ ^(4). We can now reproduce 
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block-wise the conditions given by the periodic conditions, see Eq. (2), iden¬ 
tifying the endpoints 01 , 61 , 02,62 in different ways. For instance, we can 
consider the following unitaries [ 11 ]: 


/O 

1 

0 

0\ 


/O 

0 

0 

1\ 

1 

0 

0 

0 


0 

0 

1 

0 

0 

0 

0 

1 

) U 2 = 

0 

1 

0 

0 

VO 

0 

1 

0^ 


VI 

0 

0 

0/ 


The former corresponds to the topology of two circles whereas the latter 
corresponds to the topology of one circle. 



Figure 1: Identifications of the endpoints corresponding to different choices 
of the boundary conditions. Dotted lines correspond to the topology of two 
circles, Ui. Solid lines correspond to the topology of one circle U 2 ■ 


Example 5 (Free particle on the disk, VL = D). In this case the boundary 
is given by the unit circle d^l = . The unitary group at the boundary 

is therefore given by the infinite dimensional unitary group U (£^(S'^)) . To 
impose different topological structures we shall introduce a splitting of the 
boundary dQ = Ti U r2 where Ti and r2 are disjoint and diffeomorphic 
to each other (see Fig. 2). The isomorphism £^((9D) ~ £^(ri) © £^(r 2 ) 
allows us to define unitary operators adapted to this block-wise structure. 
For instance one could select 


U = 


V 0 


0 ^ , 

l£2(r2)y 


what would impose Dirichlet boundary conditions on Ti and Neumann 
boundary conditions on r 2 . 



Figure 2: The boundary of the manifold, cID = 5^, is split into two parts, 
T 1 and T 2 . 

We want to be able to describe different topologies in this situation. From 
Fig. 2 it is clear that one cannot directly use the structure given in (2) since 
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Fi and r 2 are diffeomorphic as differentiable manifolds but not necessarily 
isometric. To overcome this difficulty we can do the following. Select a 
reference manifold Tq, diffeomorphic to Ti and r 2 but otherwise arbitrary. 
We denote the diffeomorphisms by : Tj —Tq . Fix a volume element, 
dfj,o, for the manifold Fq . Now we can define the following unitary operators 
T, = C^{Ti) ^ CHTo,dfio) 

T.ch := o g,) , (4) 


where \ Ji\ is the Jacobian determinant of the diffeomorphism gi and Hi 
the proportionality factor: g*dgri = /Ujd^o > cf. [8]. Using these unitary 
operators we can define the following unitary operator on 

0 tIt2\ 

tITi 0 ) 



This operator leads to the boundary conditions 


^i(<Flri) = T2{(f\T2) , ^i(<Flri) = -T2{(P\t2) , 


which can be called generalised periodic boundary conditions. 


3 Boundary Dynamics 

As stated at the beginning of this article we will discuss dynamics of a 
time-dependent quantum mechanical system, i.e. the Hamiltonian operator 
will depend explicitly on time. This means that we will have a family of 
self-adjoint operators H{t). In a conventional time-dependent problem, one 
typically has a time-dependent perturbation modifying the Hamiltonian cor¬ 
responding to the free particle. Even more generally, it is usually assumed 
that the functional form of the Hamiltonian operator depends explicitly on 
time. In contrast, we aim to modify the boundary conditions externally and 
leave the functional form of the Hamiltonian operator unchanged. As we 
have seen in the previous section, this amounts to change the self-adjoint 
extension of the operator. Hence the Hamiltonian operator will have a fixed 
functional form but its domain will depend on time through the boundary 
conditions. More concretely, we will have a set of parameters, u G M, 
characterising different boundary conditions that will depend explicitly on 
time. The space of parameters M will be in general a subset of the set of 
all possible self-adjoint extensions of the Hamiltonian operator. 

It is time now to discuss the time-dependent Schrodinger equation 

= Himt). 
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Recall the discussion of Section 1 where the solutions of the time-independent 
Schrodinger equation where obtained in terms of a strongly continuous one- 
parameter group of unitary operators. Similarly, solutions of the time- 
dependent equation can be obtained in terms of what is known as a unitary 
propagator. 

Definition 6. A unitary propagator is a strongly continuous two parameter 
family of unitary operators U :MxM—that satisfy the properties 

U{t,s)U{s,r) = U{t,r) , U{t,t)=ln. 

Unlike the time-independent case, there is no general result like Stone’s 
theorem guaranteeing the existence of a unitary propagator for a given time- 
dependent Hamiltonian. Hence, given a family of self-adjoint operators 
G / C M} one needs to look for sufficient conditions for it to exist. 
The situation is particularly involved since the domains of the operators are 
going to change with time. We will take now some simplifying hypothesis. 

Hypothesis HI. Let T)[u) denote the domain of the self-adjoint extension 
of the operator H defined by the boundary condition determined by n G Ad. 
For every u G Ad, the spectrum of {H,'D{u)) only contains eigenvalues 
with finite degeneracy. That means that we do not consider situations with 
singular nor continuous spectrum. This implies that for each n G Ad there 
exists a complete orthonormal set of eigenvalues 

{^wlneN • 

Under these conditions a reference self-adjoint extension (H, T){uq)) with 
orthonormal base {<hg}„gpj can be hxed. We define a family of unitary 
operators Vu ■'H ^ Ti hj 

This definition characterises completely the unitary operators Vu ■ We will 
also require the following. 

Hypothesis H2. For each n G Ad the unitary operator Vu must be a 
bijection 

Vu : ^{u) 'D{uo) . 

Notice that this is not an immediate consequence of the definition of Vu ■ 
In fact, this hypothesis imposes asymptotic conditions on the sequences of 
eigenvalues of the family {H,V{u)). 

Under conditions HI and H2 the operator VuHVu is a self-adjoint oper¬ 
ator with domain V{Uq) . Let ;= Vu^ G V{uq) for any T G V{u). Now 
we can compute Schrodinger’s equation in terms of to get; 


(6) 



To obtain this equation we have used Leibniz rule. The derivatives have to 
be understood in the sense of the strong topology. In the right hand side 
vj stands for the derivative of the family of unitary operators where the 
parameters u = u{t) are assumed to be time-dependent. This derivative will 
be in general an unbounded operator and might even not exist depending 
on the particular problem. In order to guarantee the existence of a unitary 
propagator that solves the Schrodinger equation we can use the following 
particular form of a theorem due to J. Kisynski, cf. [12]. 

Theorem 7. Under the hypotheses HI, H2 and: 

Hypothesis H3. (T , <I>) is twice continuously differentiable for all T G 

n. 

Hypothesis H4. (T , VuHVu^) is continuously differentiable for all T G 7^ 
and for all T G T’(uo) . 

Hypothesis H5. {H,V>{u{t))) is semibounded from below uniformly. 

There exists a unique unitary propagator solving the time-dependent Schro- 
dinger equation defined by the time-dependent Hamiltonian H. 

From the form of Eq. ( 6 ) one can see why the regularity and differ¬ 
entiability conditions are needed. The semibounded condition is a further 
spectral condition on the family of self-adjoint operators. In fact, the condi¬ 
tion considered originally by Kisynski is more general and includes this one, 
that turns out to be enough in our case. Under certain conditions on the 
family of unitary operators U describing the boundary condition and there¬ 
fore parametrising the different self-adjoint extensions, the Laplace-Beltrami 
operator can be proved to be semibounded. We refer to [ 8 ] for a detailed 
analysis of this situation. 

Definition 8. We say that a unitary operator U (£^(917)) has spectral 
gap if one of the following conditions hold. 

• 1/ -|- I is invertible 

• { — 1} is an element of the spectrum of U and is not an accumulation 
point. 

Theorem 9 (Ibort, Lledo, Perez-Pardo, [ 8 ]). Let U be a unitary operator 
with speetral gap. Then the self-adjoint extension of the Laplaee-Beltrami 
operator defined by U is semibounded from below. 

The lower bound is proved to be directly related to the gap of the unitary 
operator at the boundary, that is, the distance between { — 1} and the closest 
element of the spectrum to {—1} . In order to have a uniform semibound 
for the family of self-adjoint extensions it is enough to have a uniform gap 
for the family of unitary operators. 
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Example 10 (Laplacian with quasi-periodic boundary conditions). Let us 
consider the Laplace-Beltrami operator H = on the interval O = 

[0, 27r] . We consider the one-parameter family of self-adjoint extensions 
given by quasi-periodic boundary conditions defined in Ex. 3 by the unitary 
operator (3). This family of operators satisfies condition H5 since all the 
elements have spectral gap. This is so because the spectrum of every element 
of the family is cr{U) = {—1,1}. The eigenvalues and eigenvectors can be 
computed explicitly and one gets 

E” = (n + ef , , n G Z . 

\/27r 

The family of unitary operators becomes in this case 

:= . 

That is, a family of unitary multiplication operators acting on £^([0,27r]). 
It is clear that HI and H2 are also satisfied. Further calculations show that 

Hence, the regularity conditions H3 and H4 hold if the first and second 
derivatives of the function e : M —[0, 27r] are bounded. Substituting all the 
expressions into the time-dependent Schrodinger equation (6) it is straight¬ 
forward to get 

-e]^^ + MeT , T G P(e = 0) . 

This equation corresponds physically to the quantum Faraday law. It rep¬ 
resents the motion of a particle spinning in a planar wire surrounding a 
magnetic flux of intensity e that changes with time. The derivative of this 
magnetic flux is proportional to the electric potential that the particle ex¬ 
periences. 

4 Dynamical topology change 

In Section 2 we have shown how we can use boundary conditions to de¬ 
fine self-adjoint extensions of the Laplace-Beltrami operator, whose domains 
present different topological structures. In Section 3 we have addressed the 
problem of existence of solutions of the time-dependent Schrodinger equation 
when the time dependence is located strictly at the boundary conditions. 
We also provide a set of sufficient conditions (HI to H5) for solutions to 
exist. Now we will study if we are able to paste everything together to 
achieve a dynamical variation of the topology. Let us work a further exam¬ 
ple and check whether or not these conditions are met. We shall consider 
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Figure 3: Disk with boundary conditions split into four diffeomorphic re¬ 
gions. Identification along the solid lines correspond to the topology of a 
sphere. Identification along the dotted lines corresponds to the topology of 
a torus. 

the same construction introduced in Ex. 5. This time we split the boundary 
into four different pieces, /i, I 2 -, I 3 and (see Fig. 3). We can use the iso¬ 
morphism (4) constructed in Ex. 5 to identify the pieces in different ways. 
The unitary operators defining the boundary conditions will have block-wise 
structure with blocks given by (5). We consider two limiting cases for the 
time-dependent dynamics: 

• Topology of a sphere. 

• Topology of a torus. 

We need to select a path in the space of self-adjoint extensions that con¬ 
nects these two cases. Regarding Hyp. H5 it is easy to check, cf. [13, 
Corollary 3.5], that when the spectrum of the unitary operator is of the 
form cr{U) = { — 1,1} the spectrum of the Laplace-Beltrami operator is 
non-negative. Moreover, in this situation one can use spectral bracketing 
techniques, cf. [14], to show that the n — th eigenvalue of any extension 
parametrised by an unitary operator with spectrum <j{U) = {—1,1} satis¬ 
fies 

< A„(n) < A^ , 

where and A^ are the n-th eigenvalues of the Laplace-Beltrami opera¬ 
tor with Neumann boundary conditions and Dirichlet boundary conditions 
respectively. This shows that HI and H2 hold. Now we need to address the 
regularity conditions. A sufficient condition for H3 and H4 to hold would be 
that the eigenvalues and eigenfunctions were smooth functions with respect 
to t and that the norm of the first and second derivatives could be bounded 
by a constant that does not depend on n. This is exactly the situation in 
Ex. 10. Whether or not this condition is achieved will depend on the path 
chosen. The existence of such a path is part of actual research and we leave 
the analysis for future work. 
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